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Abstract 

It is widely believed that hadronisation leads to 1/Q corrections to e^e~ event 
shapes. We show that there are further corrections, proportional to {\d.Q)^/Q with 
A = 4Ca/Po — 1-6, associated with hadron masses and whose relative normalisations 
can be calculated from one observable to another. At today's energies these extra 
corrections can be of the same order of magnitude as 'traditional' 1/Q corrections. 
They fall into two classes: universal and non-universal. The latter can be eliminated 
by suitable redefinitions of the observables. 
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1 Introduction 



For many non-perturbative aspects of QCD, there exist well defined methods of study such as 
lattice calculations or chiral perturbation theory. But one of the most fundamental aspects of 
non-perturbative QCD at high energy colliders, namely hadronisation, the relation between 
what is calculated perturbatively (parton level) and what is actually measured (hadron level) 
remains well beyond the reach of such methods. At best there are complex models which 
appear to reproduce much of the experimental data, but they leave something to be desired 
in terms of a fundamental understanding of what is involved in hadronisation. 

Over the past few years there has been considerable theoretical p|^T8| and experimental 
r9|-|2^ interest in the study of hadronisation contributions to e'^e" and DIS event shapes. 



These variables provide a convenient 'laboratory' for such studies because the hadronisation 
effects are responsible for a significant fraction of the observable, making it feasible to carry 
out quantitative tests of the theoretical predictions. 

A number of the current theoretical approaches are based on the idea of extending the 
reach of perturbation theory to very low scales — thus hadronisation effects are argued 
to be related to the infrared behaviour of the coupling, or equivalently to the high-order 
behaviour of the perturbative series. All the methods predict that the hadronisation (or 
power) corrections to event shapes should scale as 1/Q, where Q is the hard scale of the 
process. They also aim to predict the relative normalisation of the corrections from one 
observable to another. 

As a result of their perturbative origins, the theoretical calculations are usually carried 
out with the assumption that all particles are massless. But in practice the observed hadrons 
do have masses. So in this paper we examine how the treatment of hadron masses modifies 
one's expectations about power corrections. 

It is perhaps worth illustrating how masses can affect our observables with a simple 
example. Let us consider the event-shape variable p, i.e. the squared invariant jet mass, 
normalised to (for brevity we often refer to it simply as the jet mass). Initially one might 
consider a Born event consisting of a pair of back-to-back particles, each of mass m. The jet 
mass has the value p = m?/Q'^. If this were the end of the story then we would argue that 
particle masses can be neglected since they give a correction. 

But one also needs to consider events containing soft particles. Here the situation is 
qualitatively different: the jet mass is [E'^ — P^)IQ'^i where E and p are the total jet energy 
and 3-momentum respectively. For a jet aligned along the particle with mass m^, 

2;-momentum p^i and energy Ei contributes an amount (£', —Pzi)IQ to the jet mass and this 
difference includes a piece of order m'f / {2EiQ) . For a soft particle with ~ A this translates 
to a 1/Q contribution. So parametrically at least, mass-related effects from soft particles are 
of the same order as traditional power corrections, and should not be neglected. 

This is not the end of the story, for the simple reason that an event generally con- 
tains many soft particles: there will be a multiplicity-related enhancement of the m?/QEi 
correction. The sum, over all particles, of l/Ei is proportional to the (—1)**^ moment of 
particle energy fractions. It has been known for a long time that positive moments and 
even the zeroth moment of particle energy fractions undergo logarithmic scaling violations 
with perturbatively calculable anomalous dimensions. It turns out that such analyses can 
be extended to negative moments, with the result that the sum 1/ E^ scales as (InQ)^, 
where A = ACa/ Po — 1-6, with /3o = (HC'a — 2nf)/3. This means that the formally 1/Q 
contribution from mass effects is enhanced by a factor (InQ)^. 

The jet mass is an example of a variable in which it is quite straightforward to see that 
there are mass effects; in other cases the mass dependence can arise more subtly. For a 



general variable, the ability to factorise transverse and longitudinal degrees of freedom (with 
respect to the quark-antiquark directions) is an essential element of traditional, 'perturbative' 
approaches to hadronisation corrections. We will discover that for an ensemble of massive 
particles, differences arise between the factorisations applying to the event-shapes on the one 
hand and particle production on the other — and by studying this mismatch we become 
sensitive to hadron mass effects. 

We will also see that a given event-shape variable can be defined in variety of ways 
(schemes) which are all equivalent for an ensemble of massless particles, but differ if there are 
massive particles. Two examples are a definition in terms of just 3-momenta (p-scheme), and 
a definition in terms of energies and angles (ii^-scheme). It turns out that of these schemes, 
one, the i?-scheme, is privileged, because it has the property that there is no mismatch 
between the event-shape and particle-production factorisation properties, even for massive 
ensembles. 

Does this mean that in the ii^-scheme there are no mass-dependent contributions? The 
answer is no: it just means that any mass-dependent contribution is proportional to the 
same coefficient cy as the 'traditional' non-perturbative correction. One can demonstrate 
that there really are still mass-dependent corrections by studying the difference between an 
event shape before hadrons have had the time to decay and after they have decayed. 

We shall present the different elements of our analysis as follows. In section ^ we define 
the event shapes that will be considered and then quote their soft limits, being careful to 
leave in the leading dependence on particle masses. 

Then in section ^ we review the tube model for non-perturbative effects 0,^, which is 
known to reproduce all basic predictions for universal power corrections. However in contrast 
to previous treatments we leave in mass effects and see how they lead to non-universal l/Q 
power corrections. This lays the ground for introducing alternative definitions of event shapes 
(massive, p, E and decay-scheme definitions) which are identical for massless ensembles but 
differ in the treatment of particle masses. 

Whereas the tube model is adequate for describing normal 1/Q power corrections it cannot 
address issues such as the scaling violations of energy moments which must be considered 
in order to derive the full Q-dependence of the mass-dependent corrections. For this we 
need to recall how coherent branching ||2^ affects hadron multiplicities and relate this to our 
particular problem. This is done in section ^. 

Then in section ^ we start to study some of the practical aspects of mass corrections — for 
example we examine how, numerically, they compare to other contributions. We discuss some 
of the issues related to the experimental measurement of mass effects (a fairly difficult task) 
and we compare our predictions to results from Monte Carlo event generators. There is good 



agreement with Herwig and Ariadne ||26|, but significant disagreement with Pythia ||2^ 
concerning the energy dependence of mass effects at very high energies. 

Finally in section ^ we use a Monte Carlo event-generator to correct data for event shapes 
to a variety of schemes and see how this affects the fitted values for the perturbative and 
non-perturbative parameters ag and a^. We also investigate the feasibility of carrying out 
fits with an extra parameter intended to allow a separation of hadronisation into 'traditional' 
hadronisation effects and mass-related effects. 

Our conclusions are presented in section ^ The appendices contain a summary of the 
notation used and introduced in this article (appendix ^), some calculational details related 
to the jet broadenings (appendix 0), some considerations about 1/Q^ power corrections 
associated with heavy quark decay (appendix |C[), and a Monte-Carlo study of some peculiar 
features of the heavy-jet mass (appendix 0). 



2 Event shapes 



The basic event shapes that we shall consider are the thrust T, the invariant jet mass p, the 
C-parameter and the total jet broadening Bt- For an ensemble of particles with momenta 
ki they are defined as follows: 

r = max2i^, (2.1a) 
(E^k^e{k,.nT)y 

(2.1b) 

3 Eo J^il 1^1 1 sin^ 9ij 
C^- '^'^ " " , ' , (2.1c) 



(2.id) 
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where is the angle between particles i and j, kj_ = {kQi.ki) the four- momentum, and ku 
the transverse momentum of particle i with respect to the thrust axis ny. 

We shall also discuss the heavy-jet mass ph which measures the invariant jet mass in the 
heavier of the two hemispheres (separated by the plane perpendicular to the thrust axis), 
and the wide-jet broadening Bw, which measures the jet-broadening in the wider of the two 
hemispheres. 

For the purpose of calculating power corrections we will be interested in the behaviour 
of the event shapes variables for two-jet configurations, since it is such configurations which 
are the most frequent, representing a fraction 1 — C (ctg) of all events. In the two- jet limit 
all variables tend to zero, except the thrust which tends to 1. Accordingly it will be more 
convenient to refer to r = 1 — T. 

Essentially a two-jet event consists of a pair of back-to-back hard particles (which may 
have fragmented collinearly) and a bunch of accompanying soft particles. The hard particles 
(associated with the qq pair at parton level) define the axis of the event, while the soft 
particles (gluons at parton level) give the deviation of the event shape variable from its Born 
value of zero. Variables like r, p and C are particularly simple (linear) in that, in the two-jet 
limit their value is given by the sum of independent contributions from each soft particle: 

^-^ E kt^e-\"^\, (2.2a) 

^ iesoft 

P ^ ^ E ^4 + ^? e-^' 0(^0 , (2.2b) 



iGsoft 



where f]i is the pseudorapidity of particle i with respect to the thrust {z) axis, ?7i = — In tan 6*^/2 
and Tji is its rapidity rii — \ \n f^j^^f^. For massless particles the pseudorapidity and the ra- 
pidity are identical. In these expressions we have neglected mass effects in the denominators 
(in the case of r, only after going from T to r), since on average the numerator is small, 
O (cis)) and thus modifications of O (A) to the denominator give effects of O {agA/Q), which 
we can neglect. 



Other variables like the heavy-jet mass and the broadenings are more complex. To see 
why, let us consider the case of the heavy-jet mass: when the event contains exactly one soft 
particle the heavy-hemisphere is always that in which the soft particle is present and the soft 
particle always contributes to the heavy jet mass. But if there are many soft particles then 
the hardest one (specifically, the one with the largest kue~^^^^) determines which hemisphere 
is heavy and any given softer particle contributes when it is in the heavy hemisphere i.e. only 
half the time. Thus the heavy jet mass is not equal to a linear combination of independent 
contributions from each soft particle. 

For the purposes of studying non-perturbative corrections we can however make a sim- 
plifying approximation: the 'hardest' soft particles come from perturbative emissions, while 
non-perturbative emissions will be much softer. The contributions from these softest particles 
do combine linearly so that one can write 



'^'^^^ E \/kl + m^e-''^e{v^) + ..., (2.2d) 



iSsoftcst 



where we have taken the hard hemisphere as being that with rj > and the dots indicate 
the contribution from harder particles. For both jet broadenings analogous arguments apply 



12], and one has 



softest 



^ J2 kuQim) + ..., (2.2f) 

is softest 



where again we have taken the wide hemisphere as being that with rj > 0. 



3 Power corrections: the tube model 

A variety of approaches exist for the study of power corrections in event shapes The 



simplest, which reproduces the results of the more sophisticated methods jToHl^, is the tube, 
or longitudinal phase-space model P,p^. 

The principle behind the tube model (the ideas of are analogous, though presented in 
a more formal language) is as follows. Soft hadrons {i.e. the hadrons responsible for non- 
perturbative corrections) are generated from the qq pair of sources; since both sources are fast 
moving (in opposite directions) a moderate boost along the qq direction still leaves us with two 
fast-moving particles, and so does not change the structure of the low-transverse momentum 
fields at central rapidities. As a result soft particle production at central rapidities must be 
boost-independent. Since a boost just corresponds to a shift in rapidity, this is equivalent to 
saying that non-perturbative (and in general soft) particle production is rapidity-independent, 
at least for rapidities rj <^\nQ/kt, beyond which one becomes sensitive to the finite energy of 
the source. The tube model makes no statement about the transverse-momentum distribution 
of soft particles, so we just write the distribution of non-perturbatively produced[| hadrons 
of type h as being: 



"'^ 'Non-perturbatively produced' is a rather awkward term — in section ^ we examine in detail what we 
mean. 



with (phikt) some a priori unknown function. There is no Q-dependence in (phikt) since the 
fields generated by a source {q or q) moving close to the speed of light are independent of 
the energy of the source. 



3.1 Massless case 



Most power correction analyses work within the approximation that particles are massless.Q 
As a result in eqs. (^]2|) all explicit mass- dependence disappears, and one can replace fj with 
rj. This leaves the expressions for the event-shape variable V in a factorised form, whereby 
all rapidity dependence (which differs from one observable to another) can be separated from 
the transverse momentum dependence {kt/Q in all observables) : 



iesoft ^ 



(3.2) 



For example for the the thrust we have friv) =6 The non-perturbative contribution to 
the mean value of the event-shape is then given by 



(Vnp) 



dkt kt 



where we have defined a non-perturbative parameter 

/ dkt^(j)h{kt) 



Oo 

Q 



cv 



(3.3) 



ao 



and a calculable, variable-dependent coefficient. 



(3.4) 



Cv 



(3.5) 



The factorisation of rapidity and kt dependence is the prerequisite for universality, namely the 
fact that the power corrections for a range of observables all depend on the same, universal, 
non-perturbative quantity (ao), with a calculable coefficient (cv). The predictions for the 
coefficients are given in table |1[ 

The more complex form for the broadenings |jl5|] arises because these variables are sensitive 
to the mismatch between the thrust axis and the quark axis. Emissions are uniform in rapidity 
with respect to the latter, whereas the broadening measures kt with respect to the former. 
After considering recoils one finds that there is an effective cutoff on contributions from 
rapidities (with respect to the thrust axis) rj > InQ/pt, where pt is the transverse momentum 
of the quark with respect to the thrust axis. This leads to cb ~ {\nQ/pt) which is of the 
order of 1 / ^/a^. 



^Even the 'massive-gluon' calculations make this approximation, since they generally assume that the 
massiv e glu on decays into two massless particles. An exception is to be foimd in as discussed briefly in 
section 3.5. 
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Table 1: coefficients of 1/Q power corrections; Po = ^-Ca — l^^f and r/o — 0.13629. 



More sophisticated approaches to the problem of power corrections at first sight seem quite 
different from the tube model: they examine the high-order behaviour of the perturbation 
series, or the dependence of the observable on a dispersive gluon virtuality (see for example 
But it turns out that both of these procedures are equivalent to determining the 
dependence of the observable on infra-red properties of the coupling; since the production of 
gluons with a given kt is rapidity-independent and proportional to as(^t) we have a situation 
very similar to the tube model but with (f)h{kt) replaced with as{kt). Accordingly we can 
write a relation between oq and the quantity ao(/U/), often used in phenomenological analyses, 
defined as 



aoifii) = I —as{kt), (3.6) 
f^i 



namely |,|,|TT],|T1 



ao 4C^^.^^L,(^,)_«.(g)-/3o#fln^ + ^ + lU (3.7) 



with M ~ 1.490, K = Ca - ~|^f- The purpose of the as terms in (|3^ ) is to subtract 
out contributions that are already taken into account in the perturbative calculation of the 
mean value. 

It has become a standard procedure to carry out simultaneous fits for and ao in mean 
values of a variety of event shapes. One important test of this class of models for hadronisation 
corrections is then that the fitted values ao (and as) should be the same for all variables, i.e. 
that ao should be universal. 



3.2 Including mass effects 

From the point of view of the tube model itself nothing changes when one introduces masses 
for the hadrons — we still have a distribution of hadrons independent of rapidity and with 
some unknown dependence on kt. What does change is that we need to use the full (massive) 
expressions for the values of the event shapes. In most cases the event shapes are just 
defined in terms of the pseudorapidities (angles) and transverse momenta of the particles. 
This means that we have to keep track of the relation between rapidity and pseudorapidity, 
which is expressed in the following (equivalent) equations: 

kt sinhf] = ^/kt + iv? sinh rj , (3.8a) 
k1 cosh^ f] + = {k^ + m^) cosh^ r] . (3.8b) 
For example for the thrust and the C-parameter, we have 

(VNP-_e) = / ^ I E / {n {v. ^) ) (3.9) 

where all that has changed compared to eq. (|3.3|) is the replacement of /v(?7) with fvijf)- 
However the fact that f/ is a function of both rj and m^/kf means that eq. ( p.9|) cannot 
be factorised into two independent pieces. Hence universality (a direct consequence of the 
factorisation) is broken. 

For a general variable V it will be convenient to write the resulting non-universal mass- 
dependent piece of the power correction as 

(5™V) = j J dv5fv(^V.^) , (3.10) 



where for example, for the thrust 



m 



5fr{ V, ^2 



frin) - friv) = e-l^l - e-l'^l 



(3.11) 



We then proceed in a manner analogous to that in the massless case. We define a Scy, 



6cv { ^ 



(3.12) 



which unlike cy depends on kf (a consequence of the non-factorisability) . Then the mass- 
dependent correction to the mean value of the event shape is 



/A i;\ f dkt . kt (ml 

(^.»v>=i:^y ^*<^v(^ 



(3.13) 



To understand the properties of the mass-dependent corrections we need to study the 5cv's: 

(3.14) 



5cr 



m 



dvfviv) 



dr] 



where 



drj 

df] 



cosh?7 



cosh^ 1] + m? /kl 



(3.15) 



For the event-shapes in the first line of eq. ( |3.14| ) the integrals have been rewritten with a 
change of variable in the first term as this simplifies their subsequent evaluation. 
The exact forms for the 5cv are 



m 



^^e^ 



m 

2 

m 



37r . 



(3.16a) 

(3.16b) 

(3.16c) 
(3.16d) 



where we have introduced the shorthand = 1 -|- m? /k^^ and E and K are the complete 
elliptic integrals defined as follows, 



E(x) 
K(x) 



7r/2 
7r/2 
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(3.17a) 
(3.17b) 




Figure 1: The dependence of cy + 6cy on particle masses, shown for the four distinct groups 
of variables. 
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Table 2: The coefficients of the high-fcj behaviour of the mass-dependent power correction, 
7v. In the case of the broadenings the cy depend on as{Q), so the ratio 7v/cv is shown for 
Q = Mz. 

To see how the Scy will affect our power correction we compare them to the universal cy 
power contribution. Figure |l| shows [cv + 5cv(m^/fct )]/cv for a range of variables as a function 
of kt/m. One immediately sees that for small kt the non-perturbative correction to the jet 
mass will be enhanced, while the NP correction to the other variables will be suppressed. 
Furthermore the enhancement for the jet mass is much larger than the suppression for the 
other variables (which are fairly similar to one another). 

To study the question more quantitatively we observe that for large kt, the 5cv scale as 
m'^/kl 




(3.18) 



with the 7v given in table |^. This means that the integral in eq. (p.l3|) is dominated by low 
momenta for all reasonable forms of the distribution of particles (phikt), and hence just gives 
a number. Therefore (5m V) is proportional to 1/Q, i.e. formally of the same order as the 
universal power correction. 

The question of the quantitative relationship between the sizes of the mass- correct ions 
in the different observables is more delicate because it depends on the region of kf which 
dominates the integral in eq. ( |3.13| ). If (pikt) is such that moderate fc^'s dominate [i.e. where 
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Figure 2: The observable-dependent piece of the integrand for the non-universal mass- 
dependent power correction, eq. (|3.13|) , normalised to 7v. The curve labelled 'decay' will 
be discussed in section E^. 



5cv is equal to 7vm^/A;^) then we can expect the following relation to hold: 

(5„V)(X7v. (3.19) 



If (f){kt) is such that smaller kts dominate the integral (p.l3|) , then formally we can make no 



such statement. Nevertheless by examining figure 0, which shows kt/mdcy/'fy as a function of 
kt/m, one observes that for r, C and the broadenings the shape of the 5c-\;{rn? /kf) functions 
are very similar. This means that regardless of the form of 0(fct) we will still observe the 
property that 

^l^^MA. (3.20) 
7Vi 7V2 

For the jet masses on the other hand, bc\i{vn? j k^) has a different shape and being considerably 
larger at small momenta, so that we can expect the following relation to hold: 



Ip 7r 

The conclusion of this section is that mass effects introduce extra 1/Q power corrections, 
which break the simple 'universal' picture of power corrections that is obtained in the mass- 
less case. For most variables we expect a negative correction, whose magnitude is roughly 
proportional to 7v (which for these variables is roughly — cv/4). For the jet masses we expect 
positive corrections whose magnitude is larger than what would be expected from a simple 
proportionality to 7v (which itself is cv/2) . 



3.3 Alternative schemes 



So far we have used the event-shape definitions given in eqs. From the point of view 

of the perturbative QCD calculation however we could have chosen any number of related 
definitions with the same massless limit and we would have obtained the same perturbative 
(and universal non-perturbative) predictions. Here we discuss two particular examples of 
such modifications. 



The p-scheme: The difference between the jet masses and the other variables occurs be- 
cause the jet masses are the only variables to be sensitive to the difference between the 
energies and 3-momenta of the particles. However one could equally well consider a second 
pair of variables, identical to the jet masses except that they are defined only in terms of the 
particle 3-momenta {i.e. in the definition, each occurrence of particle energy is replaced by 
the modulus of the corresponding 3- momentum). We will refer to these as the jet masses in 
the p-scheme (whereas we will refer to the default definitions as the massive scheme). 

As pointed out at the beginning of the section, from the point of view of the perturbative 
and universal non-perturbative calculations, which ignore particle masses, such a variable 
would have identical properties to the original jet mass. However, for an event consisting of 
soft massive particles its value would be 

Pp^^Yl ^*^e-l^^le(r^,), (3.22) 

^ iesoft 



rather than eq. ( |2.2b| ). Noting the similarity between this expression and eq. ( ^.2a| ), one 
obtains that 

i.e. relative to the universal power correction the mass- dependent piece is identical in the two 
cases. The use of the p-scheme makes no difference for the variables other than jet masses, 
since they are all already defined purely in terms of the 3-momenta. 

So in the p-scheme all variables should have a mass- dependent correction which is roughly 
proportional to 7v, which itself is roughly proportional to cy. Therefore universality should 
(more or less) appear to remain intact. 



The ^^-scheme: Another definition of event-shapes which is identical at the perturbative 
and universal non-perturbative level is one defined purely in terms of the energies and di- 
rections of particles, i.e. where all 3-momenta are substituted with momenta in the same 
direction but whose modulus is equal to the energy. We call this the i?-scheme. We note 
that similar definitions, in terms of energy fiow, have been suggested in the past by various 
authors on the grounds that they are closer to what is measured in a experimental 



calorimeter and that they may also allow event shapes to be expressed in terms of correlation 
function of fields. 

The expression for in the soft limit is 



1 y cosh^T^i + mf//c2 



Q ^ cosh?7j 
^ iesoft '* 



where the extra factor compared to Pp is the ratio of the energy to the 3-momentum. The 
expression for Jcp^ is then 



However, noting the form of drj/di], eq. ( 3.15| ), one sees that this is identically zero. A similar 



phenomenon occurs for the other variables, i.e. in general we have 



2 

m 



^cv. [-^J ^ • (3-26) 

In other words in the i?-scheme there is no non-universal mass-dependent 1/Q power correc- 
tion. So if one wants to study the universality of 1/Q power corrections, the best to way to 
do it is to measure all variables in the i?-scheme.0 

The p and E'-schemes are of also of interest because in principle one can measure the 
difference between a given observable in two different schemes. For example if one measures 
the difference between r in the p and E'-schemes, one expects this to be equal to {SmT). This 
can be done for various observables, after which one can verify the relations ( p.20| , p. 211) . 



3.4 Hadron decay 

Quite often hadron-level measurements are performed on particles which are unstable, though 
long-lived compared to their time of flight across the detector. The definition of the observable 
does not specify at what stage of the hadron decay chain we should make our measurement, 
so we are actually free to make it at any stage we like — as long as we specify the stage. 

This leads us to wonder about dependence of the observable on the particular hadron level 
that is chosen. It is possible quite generally to argue that redefining the hadron level should 
not affect the universality pattern. Suppose one starts off at a stage consisting of short- 
lived hadronic resonances. The boost-invariant nature of the mechanisms of hadronisation 
implies that these hadronic resonances will have been produced with a rapidity-independent 
distribution. Hadron decay is also a boost invariant process, so the decay products of the 
resonances will also be distributed in a rapidity-invariant manner. But the mean transverse 
momentum may well change because the decay of a hadron liberates energy, some of which 
may enter into the transverse degrees of freedom. Therefore the power correction to a given 
observable will increase as a result of hadron decay, but the increase will be proportional to 
Cv (as long as the observable is measured in the i?-scheme). 

Though it is perhaps more natural to treat hadron decay as a change in the distributions 
(j)h{kt), if one wants quantitative predictions it turns out to be more convenient to discuss it 
in terms of a dcyiw? /kl) contribution. Since our decay process is rapidity independent we 
can write the 5c\!^^^^^ for an arbitrary variable as 

-^J =CvXdecay (^-^J (3.27) 

^We note though one small defect of the i?-scheme, namely that the rescaled 3-momenta do not necessarily 
add up to exactly zero. Experimentally this is in any case quite common due to measurement errors, and 
so is not necessarily a major defect. However to preserve various desirable properties of the event-shape 
definitions, in the £^-scheme we choose to boost the event (by a small amount, of order A/Q) so as to place it 
in the centre-of-mass frame. One might worry that the boost itself might alter the value of the event shape, 
but it can be shown that the effect of the boost is to modify the event shape by a relative amount A/Q, so 
that the effect on the mean value is of order as{Q)A/Q, i.e. formally negligible. 



where J'^decay, which is independent of the variable, is defined by 



X, 



decay 



rrf \ _/ 1_ ^ 



(3.28) 



G children 



is the mean relative change in total transverse momentum as a result of the decay of a parent 
with transverse momentum kt and mass m. 

In general X^ecay is a fairly complicated function, because actual hadronic decays can 
involve 3-body final states where one or more of the decay products is massive. For the 
purposes of the studies in this paper we instead introduce a decay-scheme, in which all 
hadrons are artificially decayed to a pair of massless particles. At first sight this seems a 
little arbitrary, but there are two reasons why it is nevertheless of interest. Firstly, it is fairly 
straightforward to apply this scheme to a given ensemble of particles: one simply takes each 
hadron and in its centre of mass frame decays it to two particles moving in opposite directions 
along a randomly chosen axis. The second reason is phenomenological: as we shall discover in 
section p.3| if one applies the procedure to different hadron levels (say the normal hadron level 
and some 'resonance' level, earlier on in the decay chain) one finds that the decay scheme 
results in the two cases are very similar — in other words, decay-scheme results are almost 
independent of the particular hadron level from which one starts. 

In this 2-body decay scheme, Xdecay is given by the following expression 
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(3.29) 



for which we have yet to find a closed form. Its behaviour for kt ^ m is such that for a 
general variable the expansion coefficient defined in eq. (|3.18|) is 



7Vdecay " T^V 



1 



(3.30) 



The shape of the function 5cv^^^^ is similar to that of 5cp, as can be seen from figure 



3.5 Relation to massive-gluon calculations 

Many of the traditional power-correction calculations are based on the dispersive, or massive- 
gluon approach. Often however, as we have already mentioned, there is an implicit or even 
explicit [ITOHl^ assumption that the gluon decays into massless particles, leading to the 
statement of universality. 

One exception is the calculation in |P which considered the thrust using the full kinematics 
of the undecayed massive gluon (both in the numerator, where gluon decay often makes no 
difference, and in the denominator where it does make a difference). They obtained the 
result that the coefficient for the thrust should be AG rather than 2, where G ^ 0.916 is 
Catalan's constant. This seems quite strange since we have argued that a proper treatment 
of particle masses should lower the value of the thrust rather than increase it, whereas the 
analysis of suggests that the power correction increases. However the situation is subtle 



because the massive-gluon approach takes the power correction as being proportional to the 
non-analyticity in the gluon mass (or virtuahty) after integration over the whole of phase 
space, and there is a non-trivial relation between the effect of the gluon mass on the value of 
r for a given event and the non-analyticity. 



4 QCD-based analysis 

When discussing universal power corrections in section |^ we introduced the function (f)h{kt), 
the one-particle inclusive distribution for the non-perturbative production of hadrons h with 
transverse momentum kt. Within the tube model we vaguely know what we mean (the 
hadronisation associated with the low momentum fields from the q, q sources) , but in QCD 
it is quite ambiguous: after all, all hadrons are produced non-perturbatively! 

Strictly what we are interested in is the difference between our observable at the hadron 
level and the value calculated from a given order of perturbation theory. In our particular 
case the variables are sensitive to the mean transverse momentum (at a given rapidity). A 
proper definition of the difference in mean transverse momentum between parton and hadron 
levels is 

ao= [ (j2Mkuv) - J2 ) (4-1) 

^ * \ h p=q,g J 

where ^h{kt, v) is the distribution of hadrons h at a given transverse momentum and rapidity, 
and ^q^g{kt, rf) is the perturbative distribution of partons. So the quantity (phikt) should really 
be understood as being defined as follows: 

Y,4>h{kt) = J2MkuV) - Yl ^p(k,,v). (4.2) 

h h P=q,9 



In eq. ^]T| the integrals over ^h{kt,v) ^q,g{kt,v) separately would have values of the 
order of ase~'^'(5, because a fraction as of the time there can be a hard particle. But 
typical pictures of hadronisation state that the difference should be dominated by particles 
'produced at low transverse momenta,' as opposed to particles which come from the collinear 
fragmentation of a hard parton, since in the latter case the sum of kts of the hadrons should 
on average be equal to the k^ of the original parton and there will be no contribution to the 
difference ( |4.1| ). This ensures that the integral of (ph is dominated by low k^s and that it 
is roughly rapidity and Q-independent. 

But when working out mass- dependent effects, ^q^g does not contribute at all since in the 
perturbative calculation one has massless quarks and gluons (we do not consider the case of 
calculations with massive quarks) . So in the expression for the mass-dependent non-universal 
power correction, eq. ( p.9|) , we shall replace (phih) with ^hihyV)- 



4.1 Spectrum of hadrons 

To study mass effects in detail it is necessary to have some understanding about The 
simplest approach that is currently available is based on local parton-hadron duality (LPHD) 
[^^, namely the idea that on average there is a correspondence between the production of 
partons and the production of hadrons. One can then calculate the distribution of partons 
and expect the distribution of hadrons to be very similar. Using this idea the distribution of 



\ow-kt hadrons has been calculated as a function of kt and rapidity in 
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<l>{kt,r],\nQ A) ~ — ^ , , , , + In , . In ^ ^ . +■■■ 4.3 
/3o Vnkt/A po InQo/A Infct/A y 

where A is the QCD scale in some arbitrary scheme, Qo is an unknown cutoff below which 
parton branching stops. We have explicitly added InQ as an argument of $ to emphasise 
that it is now Q dependent. The rj dependence of $ is actually properly described by this 
formula only for large f]. But since most of our integrals in i] converge rapidly we will usually 
be able to ignore the 77-dependence altogether. 

The first term in the brackets in (|4.3| ) just corresponds to the radiation of a single gluon 
of transverse momentum kt from the qq pair, with intensity as{kt). This term is both rapidity 
and Q-independent. The second term comes from the coherent (or angular ordered) radiation 
of another gluon from the first gluon, with the logarithms originating from the integrations 
over the two gluon momenta. 

At lowest momenta the first term dominates. At higher momenta the second term becomes 
more important. We note that it has significant Q-dependence: this means that a piece of 
the mass-dependent correction will behave as InlnQ times some function of kt, and the fact 
that this function is enhanced at larger values of kt implies sensitivity to 6c^{m'^/k'f) in a 
region where the approximation of (5cv by 7v^«^/fc^ might be expected to work — in other 
words we expect there to be a term in the mass-dependent power correction proportional to 

Inlng/A 

as opposed to a simple 1/Q correction. 

A proper treatment requires that one take into consideration not only the first term 
to have Q dependence in ([4.3|) but also yet higher orders. This can be done via moments 



D^(lnQ) of the multiplicity distribution of particles with momentum fraction x emitted from 
a gluon at scale Q, D{x,hiQ): 

dr 

D^{Y = lnQ/A)= —x''[xD{x,Y)]. (4.4) 
Jo X 

The corresponding moment for emissions from a quark is just Cf/CaD^{Y). We will actu- 
ally be interested in the moments of the multiplicity distribution at fixed 77 (which we take 
positive), ^uj{v), 



/Hk f k \ ^ 

which is given in terms of by 



(4.5) 



(4.6) 



Y=\nQ/K-r) 



The full multiplicity moment satisfies the following equation, embodying coherence at 
double logarithmic accuracy (DLA) as discussed for example in 

(4.7) 



By differentiating both sides this can be written as a second order differential equation, for 
which an approximate (DLA) solution is 



D^(r)^D^(Fo)exp 

with 



^l/7r^(«s(l/)) 



(4.8) 



Hence the moment of the distribution at fixed rapidity is given by 

^^{r^MQ/^) = ^^D^{\nQ/^-r^)e-^^. (4.10) 

Eq. ( [4.7| ) and its solution eq. ( |4.9| ) are usually derived for the region around = 0. There 
they are known to give the correct leading term of 7(^, proportional to -s/al. The first set of 
subleading corrections in this region are also known and can be obtained within the modified 
leading log approximation (MLLA) which takes into account effects such as a gluon splitting 
into quarks, and the part of the Pgg splitting function which is finite at ^ — > 0. These 
corrections can be embodied into a modification of and give ||30|,|34| 



7.(«.) = 7.^^^(«s) + ^ 




u \ , n 16CAtts/2vr 

+ Po 



+ 16CAas/27r / uj^ + 16CAas/2vr 



(4.11) 



where a = llCA/S + 2ni/3Cl. Around u = MLLA effects give corrections to 7^^ of order as, 
i.e. suppressed by an amount y/a^ compared to the leading contribution; effects associated 
with the correct scale choice for as start at (9(as^^) and so do not mix with the MLLA 
corrections. 

For our applications we are actually interested in the region around u = —1 and it is not 
immediately obvious that we can apply the ~ derivation. One can envisage two sources 
of problems: firstly since the D^=_i moment is dominated by low momenta one might worry 
that its evolution is entirely non-perturbative. Secondly one may wonder whether the soft 
approximation of the Pgg splitting function, implicitly included in eq. (|4.7|) , is valid (for uj = 1 
for example it would not be). But bearing in mind that for negative u, 7a; (os) = —uj + O (as), 
we see that eq. ([4.7|) has an integrand nearly independent of y over the whole integration 
region (modulo powers of as{Y — y)), so that the integration is logarithmic. This means that 
it is dominated neither by the very soft (non-perturbative) region, nor by the region in which 
the splitting is hard, and as a consequence it is safe to write eq. (^I7|). However since we are 
not in a double logarithmic, but a single logarithmic region, we can only trust the first order 
expansion of 7^^: 

7^(as) = -u;- + O {al) , 0; < . (4.12) 

Pieces of order ag come additionally both from the MLLA corrections and from other uncal- 
culated sources such as the scale choice for as, which is beyond our control. This means that 
we are not able to go beyond leading order in our studies. 

Now that we have an expression for 7^ (as), we make the standard assumption that branch- 
ing only occurs above some scale Qq, so that D^^{Yq = InQo/A) = 1 and we arrive at the 
results 

_A 

DUY) ~ e--(^-^") (y) ^ , (4.13) 



and 



where A = ACa/Po ^ 1.565(1.714) for rif = 5(6). 

The and obtained so far correspond to expectations for numbers of gluons. The 
LPHD hypothesis suggests that for a given hadron species h we should have, 

<^UV,Y) = -Ar,,^a.e-(^-^o>'') (^^] " . (4.15) 

where Mh.uj is a unknown normahsation factor, which contains the information about the 
conversion of partons into a given hadron species h. It should depend on Fo,h = InQo.h/A in 
such a way as to ensure that the final result for is independent of Qo./i- The assumption 
of local parton-hadron duality implies that Mh,u is free of soft divergences, since these should 
all have been taken into account in the QCD treatment of gluon radiation. 



4.2 Application to power corrections 

The expression for the mass- dependent piece of the power correction is (c/. eq. ( p. 10] )) 

(^-^V) = E / ^ I / (r^, ^) ^h{K V, Y) , (4.16) 



which we can rewrite as 



(5™V) = E / ^ / ^/v,.(r/) ^hMv. Y) , (4.17) 



where we have defined (note the extra factor of kt/Q) 



Let us then expand the rapidity dependence of ^h-uj{Vy Y): 

$,,_4r7, Y) = <|.,,_40, Y)(^l-^^ + ---y (4.19) 

We see that rapidity dependent pieces are suppressed by powers of l/Y, or equivalently by 
powers of as- Since for most variables (the special case of the broadenings is discussed in 
appendix ^ the rapidity integration in eq. ([4.17|) converges rapidly, powers of r] do not lead 
to any particular enhancement and we can simply neglect the rapidity dependence of ^h-ui- 

/I 

^ Scv,^'^hMv = 0,Y), (4.20) 
ft 

where we have defined 

^-■-/^"^/-M^/f (^)"|fc.(f). (4.21) 



To understand the structure of eq. (|4.2CI| ) let us consider for now just the case of the jet mass 
(in its default, massive scheme), which has 



fc^,„._,-V*'£M£lz3i). (4.22) 



Q J 4v/^ 
Eq. then becomes 

^ Q J 2m A^-K \Qo,hJ \\nQo^h/^J 

where the integration contour passes between uj = and uj = 1. For sufficiently large Q the 
integrand has a saddle point close to = 1 and accordingly we consider its behaviour in that 
region :F| 



1-^ Wo, J \\nQo,h/A) 



A{2-ui) 



(4.24) 



Such a form holds in general for all the variables. If we use it to evaluate the saddle-point 
integral we obtain the following result 

^S„V)^,.J:^f,^..-I^(-^Y (4.25) 

^ QQo,h VlnQo,h/Ay 

When rrih < Qo,h, we can quite easily study the corrections to this result since the contour 
in ( [i.2^j| ) can be closed to the right, and the integral is equal to the sum of residues at 
a; = 1, 3, 5, ... . The first residue just gives our answer ( |4.25|) . The relative magnitude of the 
contribution from higher residues depends critically on the normalisation of the poles of (5cv,a; 
at u; = 3, 5, . . . (which can be worked out) and on the value of mh/Qo,h (which is unknown) 
and so cannot be determined a priori. However the energy dependence of these higher residues 
is much weaker: for example the residue at a; = 3 goes as (InQ/A)^/^. We are of course 
assuming that Mh-ui has no relevant non-analytic structure. This cannot be guaranteed, and 
for example if the distribution of hadrons goes as dk"^ / {k"^ +rn?) then we expect Mh-w to have 
a pole at u; = 2. This would lead to corrections to our results proportional to (InQ/A)^/^. 
So we expect that the mass-dependent power correction should go as 

(,„V)^^ln^f.0(A,„^/3|) (4.2.) 

where y4/3 in the second term could potentially be A/2 and where e represents the unknown 
(but formally universal) factors in eq. ( [4.25|) 



< in--^ % . (4.27) 



These results give us two distinct predictions. Firstly mass-dependent power corrections 
should have a leading piece which goes as (In"^ Q/ A)/Q, where A ~ 1.565 for nf = 5. Secondly 
the normalisation of this leading Q dependence should be predictable for all variables to 
within a new universal constant factor e which is intrinsically non-perturbative. There can 
be additional corrections which are beyond our control, but their scaling should be closer to 
that of a pure 1/Q term, and therefore at very high energies they will be subleading. 
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assuming Mh,-u to be free of non-analyticity around lo — 1. 



4.3 Absolute predictions 



Within the tube model (and all renormalon based analyses) we had a prediction that the 
leading hadronisation correction should scale as 1/Q. Yet from the arguments so far in this 
section we can see that, even in the i?-scheme, the differences between two different hadronic 
levels (related by the decay of some hadron species) will involve a term of order ihi^ Q)/Q. 
Therefore for an arbitrary hadron level the total hadronisation correction will also have a 
piece of order (In^ g)/g. 

One may well ask whether there exists a hadron level free of (In"^ Q)/Q corrections. For 
example if one reconstructs the various hadronic decays so as to arrive at the level of the 'first 
hadronic resonances created' then one is free of the corrections associated with hadron decay. 
But it is difficult to define what is meant by the first hadronic resonances, since one doesn't 
know how far 'back' in the decay chain to go, especially when one reaches resonances whose 
width is of the same order as their mass: at this stage resonance decay and the hadronisation 
process become intertwined. 

There are even reasons to believe that hadronisation itself could lead to contributions 
of order (In^ Q)/g, as is illustrated by the following simplistic argument: in the same way 
that the decay of a massive hadron (mass m, energy E) liberates a certain amount of en- 
ergy, roughly of order th? /2E, in order to produce a massive hadron one needs to supply 
that amount of energy. The reshuffling of momenta associated with 'supplying this energy,' 
may well affect the mean transverse momentum per unit rapidity. After summing over all 
hadrons, this implies a contribution to the mean transverse momentum proportional to the 
(— 1) moment of the energy, i.e. to (ln^Q)/Q. Since the process at play should be rapidity 
independent, we expect that for a particular event-shape variable, V, the correction will be 
proportional to cy as was the case for corrections due to hadron decay (and with the same 
proviso concerning the broadenings) . We point out that the change in transverse momentum 
associated with the hadronisation could well be negative, if the energy that is 'reshuffled to 
produce the masses' comes from transverse degrees of freedom. 

So for a given hadronic level what we can say about the (ln'^Q)/Q part of the hadro- 
nisation corrections, is that for a variable V in a scheme S it has the formQ 

,ln^Q/A , 
(cv/x, + 7vsQ) g , (4.28) 



where is the same as e defined in section but now specific to our hadronic level 
the scale /i£ relates to the (mass-dependent) change in mean transverse momentum per unit 
rapidity coming from the hadronisation and subsequent decays of resonances to our hadronic 
level t 

^' \ h P=1,9 / 

We can make one further statement: there is phenomeno logical evidence that our decay 
scheme gives a reasonable approximation to actual hadronic decays, or more specifically that 
regardless of the initial hadronic level, the decay-scheme results are almost identical (c/. 



section f).3\ ). Accordingly, for any pair of hadronic levels i and i' we expect the following 



relation to hold 

4/i£ + e£ ~ 4/i£/ + ee , (4.30) 

where for variables other than the broadening we have exploited the fact that 4:'jy^^^^^ = cy. 
^In the case of the broadenings the fuU form actually has 47^^^^^ rather than c^^ii- 



4.4 Infrared and coUinear saftey? 

All the event-shapes considered in this paper are generally considered to be infrared and 
coUinear (IRC) safe. Yet above we have argued that they are sensitive to hadron multiplicities 
which are inherently IRC unsafe. How can these two statements be reconciled? 

Event-shapes are perturbatively IRC safe because they are linear in particle momenta; so 
if a parton of energy £'0 splits coUinearly into two partons of energies Ei and E2 then the 
value of the event shape is unchanged, V Ei/Q + E2/Q = Eq/Q. 

Mass effects behave differently because, simply kinematically, they are proportional to 
irP'/EQ. They are not usually to be seen in perturbative calculations because partons are 
considered to be massless. But hadrons are massive and since energy is now in the denomi- 
nator, the relation 1/Eq^ I/E1 + I/E2 means that mass effects appear to be IRC sensitive. 
There are also situations where one would expect multiplicity enhancements, similar to those 
discussed here, in purely perturbative calculations. With massive b quarks for example, in 
the difference between p and E schemes one would see an mh/Q suppressed contribution. 
This would be sensitive to the multiplicity of slow large-angle b quarks, which at high orders 
is be enhanced by infrared and coUinear logarithms of rrib/Q. 

But mass effects do not turn an IRC safe observable into an IRC unsafe one. They are 
always suppressed by powers of m/Q, so even if this factor is enhanced by logarithms of 
m/Q, in the limit Q 00 the net contribution still goes to zero, as required for IRC safety. 

5 Comparison to Monte Carlo 

It would be nice to test our predictions of mass effects against data, for example by looking 
at the differences between measurements of the same variable in two different schemes. The 
experimental difficulties are significant though. 

To calculate event shape observables in an experiment, four-momenta have to be recon- 
structed from the tracking and calorimetric data. As a simultaneous measurement of p and E 
is far too imprecise to constrain the mass, one usually explicitly assigns a mass to each of the 
reconstructed particles. The mass assignment is based on the signature in the detector. In 
general it allows the separation of neutral from charged and electromagnetic from hadronic 
particles, but the separation of different hadrons (vr, K, p or n) is experimentally much more 
difficult and usually left to specialised tagging algorithms. Often the pion mass is chosen to 
be assigned to charged hadrons, as pions are the most common charged hadrons. For neutral 
particles, in principle calorimetry could distinguish between electromagnetic and hadronic 
particles, but in practice this is difficult and all neutrals are assigned zero mass. 

In any case the effect of misassignment needs to be corrected using Monte Carlo sim- 
ulation. This means that 'measurements' of differences between schemes depend critically 
on the extent to which the Monte Carlo simulator gives an accurate description of features 
of the data which are not measured. Of course the Monte Carlo programs have themselves 
been tuned in order to reproduce the (separate) data on the production of different hadron 
species, but this tuning does not constrain all the available degrees of freedom. As a result 
it is difficult to establish the magnitude of those systematic errors on the measurement that 
are associated with the dependence on the particular Monte Carlo model that has been used 
for calculating corrections. 



5.1 Magnitude of mass effects 



When comparing the results from Monte Carlo simulations with data one cross-check comes 
from comparing the absolute value of an observable. Table ^ shows the invariant jet mass 
calculated in different schemes from Monte Carlo as well as from data. Within the exper- 
imental errors the simulations agree with the experimental results. The difference between 
the schemes, though, is more sensitive to the choice of the Monte Carlo program. Different 
simulations deviate up to ±10% indicating a systematic uncertainty of that order. The dif- 
ferences computed from the DELPHI results confirm that the simulation used by DELPHI 
gives consistent results. 

It is also interesting to see how the shift caused by switching from default to p-scheme 
depends on the value of the observable and how the different particle species contribute. 
Figure |^, for the jet mass, was obtained using Ariadne separating the contribution from 
different particle species by applying the p-scheme only to particles of a given type and 
stacking the results. 

Because the particle masses enter quadratically into p nucleons give the biggest effect, 
despite their low multiplicity. At low values of p changing the scheme causes a low average 
shift, which first rises with p and then falls off again. This structure stems from all particle 
types. An additional peak at large values of p is due to nucleons only and probably arises 
due to peculiarities of baryon production in 4-jet events. 

The overall structure reflects influences from the numerator and the denominator in the 
definition of p. Low values of the numerator are probable only at low multiplicities which in 
turn only have small corrections. With larger values of p the average multiplicity rises and so 
does the average difference between standard and p-scheme. For even larger values the relative 
change in the numerator decreases and mass-effects in the denominator get progressively more 
important, finally over-compensating the changes in the numerator. 

The interplay between numerator and denominator results in a complex p-dependence 
of the difference between schemes and demonstrates that mass-effects have a non-trivial 
influence on the observables shape. Mean values are, however, dominated by small values of 
the numerator, where the influence of the denominator is suppressed (see section y,.2\) . 

While the necessary experimental corrections make it difficult to measure the difference 
between the schemes from data, the experiments' reliance on simulation allows us to correct 
the existing data to any desired scheme without introducing significant additional systematic 
errors. We shall use Ariadne p6| to transform existing data [|T9l,|2T|,|36| to a desired scheme 



Observable 


Her wig 6.1 
default 


Pythi 
default 


a 6.1 

tuned 


Ariadne 
default 


DELPHI 

data 


P 
Pp 

PE 


0.0363 
0.0316 
0.0326 


0.0372 
0.0326 
0.0336 


0.0371 
0.0332 
0.0341 


0.0375 
0.0330 
0.0340 


0.0370 ± 0.0005 
0.0327 ±0.0003 
0.0335 ± 0.0003 


6pp 

SpE 


0.0047 
0.0037 


0.0046 
0.0036 


0.0039 
0.0030 


0.0045 
0.0035 


0.0043 
0.0035 



Table 3: Comparison of MC with data for standard and p-scheme at 91.2 GeV. The statistical 
errors on the MC results are below 0.0001. Differences between the different models indicate 
systematic uncertainties, which for 6pp are of the order of 10%. DELPHI numbers were 
obtained by averaging the published results for heavy and light-jet masses. The tuning 
used for the tuned Pythia results is based on [p5| . 
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Figure 3: The contributions to (5p, the difference between the mean jet mass p in the default 
(massive) and p-scheme, coming from different hadron species (vr^, kaons and nucleons), 
shown function of p in the default scheme. 
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Figure 4: The jet mass p as measured in the default scheme (upper set of points) and 
corrected to the p-scheme (lower set of points). The upper curve corresponds to a fit to the 
default- scheme points, while the lower three curves correspond to a fit to the p-scheme points 
(the dotted curve is the pure as component, the dashed curve includes the O (a^) term while 
the solid curve includes a power correction of the form eq. (|3.7|)). 




whenever corresponding measurements in this scheme do not exist. Currently data exists for 



non-standard schemes only from DELPHI [21], HI and ZEUS 



In figure ^ we show data for the jet mass p as function of Q in the default and in the 
p-schemes. As p is usually not given by the experiments, it was taken as the average of the 
measured heavy and light-jet masses. The lines in figure ^ correspond to fits using the O (af ) 
perturbative prediction and a power correction of the form eq. ( p.7|) . 

The difference between the two middle (red solid and dashed) curves corresponds to the 
normal, 'universal', power correction. The difference between the upper two (red and black 
solid) curves is the mass-dependent power correction. One sees that above Q = Mz it is 
as large as the traditional 1/Q term, and it can be as much as 10% of the mean value of 
the observable (c/. also table For a general observables differences between the p and 
i?-schemes at Q = Mz are of the order of a few percent, whereas differences between decay 
and i?-schemes are between 5 and 10% of the observable. 

So while none of our previous analysis has given us any indication of the absolute size 
of mass-induced effects, comparisons in this section have shown that the absolute size of 
mass effects at Q = Mz due to hadrons can be sizable portion of the non-perturbative 
power-corrections. 

Further results of comparisons with data, transformed to various schemes, will be pre- 



sented in section 



5.2 Comparison to predictions 

Here we compare Monte Carlo results with the predictions of sections |^ and ^ There are two 
main predictions which we wish to test. Firstly that the leading energy- dependence of mass- 
dependent effects is {In^ Q)IQ] secondly that the coefficient of the leading energy dependence 
is proportional to 7v. We are also interested in examining a third, more qualitative prediction, 
namely that certain subsets of observables have similar subleading 1/Q mass effects. 

We shall study results from three Monte Carlo event generators: Herwig |^ , Pythia ||2^ 



and Ariadne |^6[. Let us first illustrate the kind of behaviour that is seen by examining three 
observables: the difference between the p (default) and ii^-schemes for the C-parameter, the 
difference between the massive (default) and ii^-schemes for the jet mass and the difference 
between the decay and ii^-schemes for the thrust. These differences (multiplied by Q) are 
shown as a function of Q in figures ^a, and ^ for Herwig and Pythia and Ariadne 
respectively.^ They have all been normalised to the appropriate 7v. 

Pure 7v/Q corrections would lead to superimposed flat hues. The fact that the lines rise 
for all three event generators is consistent with the fact that we have a correction enhanced 
at larger values of Q. But the nature of the Q dependence is not consistent between the 
different programs. For Herwig and Ariadne the second derivative is positive and roughly 
consistent with In^Q with A ~ 1.6 as predicted in eq. ( [1.26| ). On the other hand the Pythia 
results are inconsistent with such a hypothesis — we return to this problem shortly. 

Our second prediction in eq. ([4.26|) was that the Q-dependence should be proportional 
to 7v — the fact that our observables (normalised to 7v) have very similar Q dependences 
supports this hypothesis. This is true for all three event generators. 

To study these questions more systematically we fit a formula of the following formQ 

(5-V) = ^f6ln^-? + b] , (5.1) 



Q \ A 

to (^mV((5)), where e, Aefr and B are the fit parameters (we take A = 0.2 GeV). Implicit in this 
procedure is the assumption that terms with subleading energy dependence are reasonably 
well approximated by ^^BjQ (fits involving more sophisticated forms for the subleading 
terms turn out to be fairly unstable). To reduce the impact of subleading effects we only fit 
points with Q > 100 GeV. We have generated 10^ events per point. 

We expect Ags to be somewhere between the rif = 5 and nf = 6 values of A = 1.565 and 
A = 1.714. The results for A^s are shown in the left hand plot of figure for all three event 
generators, together with bands representing the predicted rif = 5 and rif = 6 values for A. 
Almost all the Pythia results have A^s ~ which is a signal of a dependence of the form 
InlnQ/A. The Herwig and Ariadne results in general lie close to the predicted value for A. 

^The plots have been generated using only events with primary down-quarks: there seem to be slight 
differences between the results coming from different light-quark species and the change in flavour composition 
of events as Q approaches Mz leads to an extra small but spurious Q-dependence. More important though 
is the removal of top-quark production: for Q > 2mt the structure of the Born level of primary top-quark 
events is very different from that primary light-quark events, because of the top decay. 

^In the case of the broadenings in the decay scheme we actually use a more complicated form, in line with 
the discussion in appendix ^: 

(<5mV) - ^ (^7v,..=ieln^°"' | -f 7v,^=ooS 
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Figure 5: Herwig, Pythia and Ariadne results for the differences between pairs of measure- 
ment schemes for three observables. In all cases the differences have been normalised to 7v. 
The lines are fits of the form eq. to the points with Q > 100 GeV. These figures have 
been generated using events with primary down-quarks; the labels {d — E), (M — E) and 
{p — E) indicate differences between the decay and E-schemes, the massive and i?-schemes 
and the p and E'-schemes respectively. 



Most of the observables lie to within ±0.1 of the expected range, the notable exceptions being 
the difference between the decay and ii^-schemes for the heavy-jet mass and the broadenings. 
The broadenings in the decay scheme are quite delicate observables because they involve 
an expansion in powers of y/a^ which may be quite slowly convergent. In the case of the 
heavy-jet mass it is not too clear what is going wrong though it may well be related to the 
non-inclusiveness of the variable (c/. appendix 0). 

It is of interest to establish whether the inconsistency between our prediction and Pythia 
is due to the nature of the hadronisation or to the parton showering. If it were the former 
we might think that we had been too naive in assuming LPHD and an absence of qualitative 
changes due to hadronisation. Two distinct arguments support the hypothesis that the prob- 
lem lies with the parton showering. One is that Ariadne uses the same string fragmentation 
routines as Pythia and so the problem must lie in the part of the physics which is treated 
differently between the two programs (the parton showering). The second argument comes 
from a direct investigation of the Pythia parton level: in Pythia, parton level gluons are 
massless so we cannot simply look at the difference between two mass-schemes. However 
mass effects are just related to the sum of the inverse energies of all the particles, so we can 
instead look directly at the behaviour of D-i{Y = InQ/A) (c/. eq. ([4.4| )) at both parton 
and hadron level (where it is best examined for individual hadron species) and check that it 
has the right energy dependence. We find that in Pythia both at parton and hadron levels, 
D_i(Y) rises too slowly with energy, while Herwig for example shows an energy dependence 
which is consistent with our prediction, both at parton and hadron level. This suggests that 
the parton showering present in Pythia might be lacking some of the dynamics associated 
with the coherent branching approach used in section ^ 

Problems with coherence should have implications also for hadron multiplicities. If we 
restrict ourselves to 'uds' primary-quark events, we find that the ratio of Herwig and Ariadne 
TT^ multiplicities is essentially independent of Q (to within 1% for Q between 100 GeV and 
30TeV). The ratio of Pythia to Ariadne vr^ multiplicities on the other hand decreases by 
about 13% over this range. We note in passing that the vr^ multiplicities from Herwig and 
Ariadne at any given Q differ by about 8%, and that if one includes heavy primary quarks 
the situation is more complicated. 

Our second prediction was that the leading Q-dependence should be proportional to 7v. 
To test this systematically we fix Agg to be equal to A = 1.565 and fit for the values of e 
and B in eq. ( p.l[ ) (using the same range of Q as before). We do this only for Herwig and 
Ariadne. The results are shown in the right hand plot of fig. ^ The mean value of e is 
about 0.023 GeV for Herwig and about 0.016 GeV for Ariadne. One should not be misled into 
thinking that these small numbers imply small effects — they get multiplied by In'^Q/A, 
which is about 17 for Q = Mz\ The range of e values for different observables is typically 
about ±10% from the central value (with the exception of the difference between the decay 
and i?-schemes for ph). Thus our two main predictions, concerning the energy dependence 
and the relative normalisation of mass- dependent effects are in remarkable agreement with 
Monte Carlo results. 

We can also examine the mass correction at a given fixed value of Q, rather than its 
Q-dependence. In the left-hand plot of figure ^ we show Q(5mV(Q))/7v for Q = Mz- The 
points seem to fall into two groups: those corresponding to differences between the p and 
E'-schemes, and those corresponding to differences between the massive and ii^-schemes (jet 
masses) and the decay and E schemes. The differences between and p and i?-schemes are all 
governed by the functions 5cv in figure |^ which have a maximum: these functions have very 
similar shapes, meaning that whatever the form of ^h{kt,y) in ( |4.16| ) the integral will be 
proportional to 7v. This statement was made earlier in the form of eq. ( |3.20| ). The fact that 
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Figure 6: The left-hand plot shows the power Aes 'measured' for a range of observables 
from Pythia, Ariadne and Herwig; the yellow bands are the expected values for nf = 5 and 
rif = 6. With Aeff fixed to its predicted nf = 5 value, the values of e fitted for the different 
observables are shown in the right-hand plot (just Ariadne and Herwig). The fits are carried 
out for Q > 100 GeV to reduce their sensitivity to subleading effects using only events with 
d primary quarks. Errors are statistical from 10^ events at each of 17 energies. 
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Figure 7: The left hand plot shows the value of the observables at a fixed value of Q = Mz, 
normalised to 7v. The right hand plot shows the Q-dependence of the various observables 
normalised to the leading-order perturbative coefficient rather than to 7v. Results are shown 
for Herwig only. 



the points for the jet masses (massive minus ii^-schemes) and the decay minus ii^-schemes are 
higher was anticipated in eq. ( p.21[ ), since after accounting for the rescahng by 7v, their 5cv's 
remain larger than those of the other variables. 

Finally, for entertainment purposes, in the right hand plot of figure |^ we show the Q- 
dependence normalised not to 7v but to Apt, the first coefficient of the perturbative expan- 
sion of the event shape. It illustrates the fact that mass effects vary significantly in size and 
sign from one observable to the next — and that the ability to predict that pattern of this 
variation is a non-trivial achievement! 



5.3 Resonance and hadron-level decay scheme results 

In section |3.4| we mentioned that a phenomenological advantage of the decay scheme is that 
regardless of the hadronic level from which we start the decay scheme results are very similar. 
To verify this we start with two hadronic levels: a normal hadronic level, as defined earlier 
and a resonance level which is 'defined' as the first level of hadrons produced in Pythia 
or Ariadne. We then look at the difference between decay-scheme event-shape values for 
these two hadronic levels compared to the difference between the £'-scheme values. Results 
obtained from Ariadne are given (in %) in table ^ they show that for most variables, in 
the decay scheme one is very insensitive to the initial hadronic level. This is not completely 
trivial since the actual decays that take one from resonance to hadron level are not just the 
artificial massless two-body decays of the decay scheme. 

For most of the variables the ratio shown in table ^ scales roughly as i.e. at these 
energies the difference between decay schemes for the two hadronic levels is dominated by a 
correction rather than a, 1/Q correction. The jet mass p and the wide-jet broadening 
seem to be more complex (note also the different sign of the correction compared to the other 
observables, and the somewhat larger value for p), but the origins of the differences have yet 
to be identified. 

A point worth noting (we will see a related point in section ^ is that at lower energies, 
the good correspondence between the two decay-scheme results holds only if heavy-quark 
decays are taken into account separately. 



5.4 Total hadronisation 

So far in this section we have examined differences between various measurement schemes and 
various hadronic levels. We observed in section [4.3| that, since the difference between any two 
hadronic levels contains terms proportional to (ln(5)^/Q, the total hadronisation corrections 
in going to an arbitrary hadronic level must also contain such terms. Additionally, hadroni- 
sation itself might introduce a contribution proportional (ln(5)^/Q, as a consequence of the 
reshuffling of momenta associated with the production of massive hadrons. We introduced 
the parameter pi to represent the normalisation of such a component for a given hadronic 
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Table 4: Percentage dependence of decay scheme results on the choice of the initial hadronic 
level ('normal hadron' or resonance level), relative to the dependence of ii^-scheme results on 
the hadronic level. The numbers are shown for Q = Mz using results from Ariadne. 
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Figure 8: The difference between liadron and parton level results for r (multiplied by Q), 
shown as a function of Q. Curves are shown for three different hadronic levels, and have 
been obtained from the Ariadne event generator with only light primary quarks. The thrust 
is defined in the ii^-scheme. 



level i, and pointed out that it could quite conceivably be negative, for example if the energy 
required to produce the hadron masses comes partially from transverse degrees of freedom. 

One may wonder what happens in the hadronisation models used in Monte Carlo event 
generators. The various curves in figure ^ show the corrections to r (multiplied by Q) in 
going from parton level to each of a variety of hadronic levels, as determined from Ariadne. 
The first level after hadronisation is the resonance level — the fact the corresponding line 
has negative slope means that fie is negative. Indeed the hadronisation corrections change 
sign at around 200 GeV! 

However we expected mass effects to have a characteristic signature, namely to contain 
a term {\nQ)^/Q with A ^ 1.6. If we carry out a fit (analogous to those carried out in 
section p.2| ) to determine the effective power for the resonance level curve then we obtain 
Aeff ~ 1.0. This could be due to sub leading mass effects, or some completely different effect 
which has yet to be considered. 

At other hadronic levels one sees a weaker Q dependence — the positive (In Q)^/Q con- 
tribution from hadron decays cancels a large part of the negative contribution from the 
hadronisation. This in itself is an interesting, and even quite natural, result: the negative 
contribution from hadronisation is of the same order as the positive contribution from the 
decay of all hadrons! 

Of course these observations may be very specific to the hadronisation model being con- 
sidered. It would have been interesting to carry out a similar exercise with the Herwig event 
generator, however there the situation is more complex because the gluons have a (large) 
mass (0.8 GeV). So in some sense, part of the hadronisation mass effects will already be im- 
plicitly included in the parton showering and a straightforward investigation of the difference 
between hadron and parton levels will show up only a part of the mass effects (and there 
will be a large ambiguity coming from the choice of scheme in which to measure the parton 
level) . 



6 Fits to data 



6.1 2-parameter fits 

Following the suggestion of it has become standard procedure in recent years to carry out 
simultaneous fits for ag and ao, with formulae of the form 

where the Ci are the perturbative coefficients for the mean value and is defined in terms 
of ao and as in eq. ( p.7|) . In order to set the scene we show in figure ^ one-a confidence-level 
contours from such fits to data [|r^, for a range of event-shape variables, all in the 



default schemes. In the absence of mass effects the universality hypothesis states that the 
values of ag and ao should be consistent for the different variables. 

Compared to the figures of this kind that one usually sees, one difference is the inclusion 
of a result for p — until now generally only ph has been studied. While data do not exist 
for p itself, there are some data on the light-jet mass pi, and from this one can calculate 
P = {Ph + Pi)/'^- What one sees is a significant inconsistency between this variable and the 
others.0 

Of course we know that we should really be carrying out the fits with additional terms 
of the form ( |4.2<j| ), so as to take into account mass-dependent corrections. Let us for the 
time ignore p£ {i.e. pretend it is conveniently zero!) and concentrate on the term involving 
7VsQ — this piece is measurement scheme- dependent. In the default schemes is positive 
for the jet masses, and negative for all the other variables (c/. table If we ignore it then 
our fit parameters for p and ph should come out larger than for the other variables. This is 
exactly what we see in figure ^a. 

So if we want to be check universality we first have to ensure that these non-universal 
mass effect are absent, i.e. choose a scheme in which the 7vg are zero, namely the E'-scheme 
(we could also use a scheme in which all the 7v are more or less proportional to cy, such as 
the j9-scheme). Accordingly in figure ^ we repeat the fits for ag and ao but using E'-scheme 
data.0 The arrows show how the best fit values have moved in going from the default to the 
ii^-schemes. 

The switch to the i?-scheme decreases the values of the jet masses, while it increases, 
by somewhat less (in accord with the opposite sign and smaller value of 7v/cv), the values 
of the other variables. We see a change in both as and ao because it is only through a 
linear combination of the as and ao Q-dependences that the fits can mimic a term of order 
(In^ Q)/Q. The 'angle' of the arrows depends on the relative amounts of (In"^ Q)/Q and plain 
1/Q in the mass-correction: if mass effects involved just 1/Q corrections then only ao would 
change. For the broadenings the situation is more complex because the 'universal' power 
correction goes as l/{y/a^Q) which is more similar to a mass effect than a pure 1/Q term, so 
there is less need for a change in as to mimic the mass effect. 

^It is perhaps ironic that th is should be the one 'standard' variable that had not been studied until now! 

^As discussed in section p?i\ t here is as yet very little data in the i?-scheme (in e+e" it exists only from 
Delphi, for the jet-masses ||2l[]), so we use Ariadne to correct data from the default scheme to the E- 
scheme. To be as close as possible to what the experiments use one might have preferred Pythia, but as we 



have seen in section 5.2 this does not reproduce the correct energy dependence for mass effects, though at 
phenomenologically relevant energies the discrepancy is quite small. All other corrections to different schemes 
and hadron levels are also done with Ariadne. It should be kept in mind that, especially at low values of Q 
there are big differences for example between Herwig and Ariadne — consider r [d — E) in figure ^ — this 
implies non-negligible uncertainties regarding the effect of the scheme changes on the best-fit values for as 
and ao- 
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Figure 9: 1-a confidence-level contours from fits to event-shape variables in a range of 
schemes, (a) fits in the default schemes (normal hadron level); (b) fits in the £'-scheme 
(normal hadron level), with arrows indicating the motion of the contour in going from the 
default to the £^-scheme; (c) fits in the i?-scheme at resonance level, with arrows indicating 
the motion of the contour from the decay-scheme, to the hadron-level ii^-scheme, to the res- 
onance ii^-scheme — here the correction to resonance level has carried out using only events 
with light primary quarks; (d) fits in the i?-scheme at resonance level where the correction to 
resonance level now includes events with heavy primary quarks as well — the arrows indicate 
the motion from the 'uds' resonance level. 



In the ii^-scheme the universahty picture changes with respect to the defauh schemes: 
whereas in the defauh scheme p was clearly inconsistent with the other variables, in the 
£^-scheme it is now very close to the thrust and the C-parameter. The heavy-jet mass on the 
other hand now seems to be the least consistent of the different variables as can be verified 
by examining the contribution from ph in a simultaneous fit to all variables. It may well 
be that the non-inclusiveness of this variable is responsible for its different behaviour, as is 
discussed in appendix 0. 

Now that we are in a scheme in which non-universal mass corrections have been eliminated 
we can turn our attention to the question of universal mass corrections, i.e. the contribution 
related to pi in eq. ( [4.281) . A first question of interest is whether we actually need to worry 
about this term at all — maybe it is sufficiently small that it can be ignored altogether. We 
know that p^ depends on the hadronisation level I. To gauge the importance of p^ we study 
three levels, each in the ii^-scheme: the decay level (actually the decay scheme of the usual 
hadron level), the usual hadron level and a 'resonance' level. The latter is taken (arbitrarily) 
to consist of the first level hadrons produced in the Pythia/ Ariadne string hadronisation 
routines. The results of 2-parameter fits to these different hadronic levels are shown in 
figure ||c: the arrows start from the decay level best fit, go to the usual hadron level best fit, 
and then to the resonance level, for which we also show the l-cr contours. 

All variables move more or less in the same direction and by the same amount — this is 
consistent with our knowledge that pi is multiplied by cy (the broadenings are more complex 
and move a bit differently). Accordingly the situation regarding universality is essentially 
unchanged (if anything, in the resonance level it is somewhat improved). However in going 
from the decay to the resonance scheme as changes by up to 0.007; ao is also sensitive to the 
hadronic level chosen and varies by up to 0.2. At Q = Mz the variation in the observables 
themselves is of the order of 10 to 15%. In other words mass effects, even in a 'universal' 
scheme, are responsible for a significant part of an observable's value and have a non-trivial 
effect on fits for as and a^. 

If we had examined the same three hadronic levels in the p-scheme we would have seen 
even larger dependence on the scheme because of the contribution from the 7Vpe£ term (recall 
that e^, which is zero in the decay level, increases as one goes towards the resonance level, 
c/. eq. ( |4.30| )). In particular the dependence of as on the hadronic level should double, a 
consequence of the relation 7Vj, — — 7Vdeca • '^^^ complete range of values for as and ao in the 
different hadronic levels and E and p schemes is summarised for our different event-shape 
variables in figure |10 . 

Coming back to figure one important point which we have yet to mention is that the 
corrections from normal hadron to the resonance level have been calculated using events with 
only light (uds) primary quarks. The corrections have then been applied to all events (in- 
cluding those with heavy primary quarks) . This is equivalent to reconstructing all resonances 
except those associated with primary heavy quarks. The reason for doing this is that in an 
event with heavy primary quarks, going to the resonance level involves a reconstruction of 
the heavy-quark hadrons. In the usual hadron level these have decayed and the invariant 
mass iTLq of the hadron contributes to the event shape at the level m?^/Q'^ (for a more detailed 
discussion, see appendix and also the discussion of heavy-quark decays in the context of 
fragmentation functions in |^), whereas in the resonance level [E or p) schemes the mass 
of a forward moving hadron has little impact on the value of the observable (other than at 
order a.^Tr?JQ'^\ If we had shown figure |c including heavy-quark events to carry out the 
correction from hadron to resonance level, the combination of the heavy-quark decay effect 
and the usual light-hadron mass effects would have made it difficult to interpret the figure. 
This is the reason why we applied a correction calculated using only events with light primary 
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Figure 10: Plots showing ctg and ao values obtained by fitting to data corrected to a variety 
of schemes and hadronisation levels. Corrections relative to the default schemes have been 
obtained using Ariadne, as discussed in the text. 
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r res. (uds) 


0.1448 ±0.0141 


0.891 ±0.237 


-0.077 ±0.042 


56.2/44 


r res. (full) 


0.1206 ±0.0095 


0.262 ±0.161 


0.012 ±0.027 


57.7/44 



Table 5: Results from a 3-parameter fit to the thrust in two different hadronic levels ('uds' 
and full resonance levels, ii^-scheme). 



quarks. 

In order to give the 'full story', in figure ^ we show the 1-a contours which come from 
including primary heavy quarks in our correction to resonance level (the arrows come from 
the 'uds' resonance level). The arrows are roughly at right angles to those in the other plots 
— this is closely linked to the fact the difference between the 'uds' and full resonance levels 
is a l/Q"^ correction {i.e. something which dies off faster than 1/Q, so that we need the 
difference between a 1/Q and an as{Q) term in order to approximate it over the given energy 
range). Despite their significant effect on the best fit for as, the heavy-quark decay changes 
the value of the observables significantly only at low Q values. At Mz for example the effect 
on the thrust is 2% (compared to 16% from the decay of all the light-resonances). 

A final point worth mentioning is the following: if one believes, as implied by Ariadne 
in figure ^ that of the various hadronic levels the decay-level is that with the smallest 
contamination from mass effects, then choosing this hadronic level for a fit to ag and ao, and 
accounting also for heavy-quark decays (by either reconstructing the heavy-quark hadron, or 
explicitly including the contribution from heavy-quark decay as calculated in appendix |Cp 
will lead us to high results for as, around 0.130. This means that there may well be room 
for the large higher-order perturbative coefficients predicted by Gardi and Grunberg in WE 
which for normal j9-scheme measurements imply a value for ag of about 0.110. 



6.2 3-parameter fits 

That fact that the event-shape values depend significantly on the particular hadronic level 
chosen implies an important contribution to the event-shape value from a term proportional 
to lie. This means that we should really be fitting for fie as well as for as and ao, using an 
equation of the form 

^^^ = ^^2^ + ^W^ "Q^ A- ^ ^ ^ 

In the last term we fix A = 0.2 GeV (tying A to the value of as{Mz), in practice the replace- 
ment of In £ by 27r/ {/3oas), makes little difference). It turns out that the inclusion of the last 
term makes the fitting procedure quite unstable, and leads to large (correlated) errors on the 
individual fit parameters as well as a strong dependence on subleading effects. Consequently 
the significance of the results is limited. 

To illustrate the point, we consider the thrust at two hadronic levels (both in the E- 
scheme): the 'uds' resonance level and the 'full' resonance level. The former is subject to the 
ml/Q^ corrections arising from heavy quark decay, whereas the latter should be free of them 
(but will still have other l/Q"^ corrections). The fit results in the two cases, table ^, give 
completely different pictures. In the 'uds' resonance level, the fit results seem inconsistent 
with our expectations for ag (if only at 2-0"). In the full resonance level the value for as is 
'as we would like', but the is a bit larger. 



There may of course be other subleading effects that we have not yet considered (for 
example higher-order perturbative corrections) which could cause further big changes. Ad- 
ditionally the event generator used to correct to a given scheme may not have the right 
behaviour at low values of Q. So the systematics are such that, at least currently, it is diffi- 
cult to extract reliable information from a 3-parameter fit: the uncertainties on as are even 
larger than those which arise by considering a range of hadronic levels in a 2-parameter ctg, ao 
fit and /z^ remains unconstrained. 

If our main aim is to determine /i^ then we can try a 2-parameter fit with fixed. 
But here again we find that the systematic uncertainties on are of the same order as the 
parameter expected size of the parameter itself, i.e. about 0.02 (c/. the values of e in figure ^). 

7 Conclusions 

The Good. In this paper we have understood many features of the contributions to event- 
shapes that are associated with hadron masses: there are two classes of contribution, both 
of which scale as {hiQ)"^ /Q. The contribution from the 'non- universal class' depends on 
the details of whether the variable is defined in terms of 3-momenta, energies and angles, 
or a mixture of the two, and we can calculate the relative magnitude of the mass-correction 
for different definition schemes. It turns out that there is a privileged scheme (i?-scheme) 
in which non- universal mass corrections are absent, because their coefficient is zero. The 
'universal' class of mass corrections gives contributions proportional to the same coefficient 
Cv that appears in calculations of traditional universal (massless) 1/Q corrections; universal 
mass corrections are present regardless of the definition of the event-shape variable, and they 
are proportional to a new non-perturbative parameter //^ which depends on the particular 
hadronic level £ {i.e. stage of the decay chain) at which we observe the event. 

In traditional power correction analyses, 1/Q power corrections are often given a quasi- 
perturbative interpretation in terms of an infrared-finite coupling, reflecting the fact that 
they are associated with the strictly perturbative concept of the renormalon divergence of 
perturbation theory. The parameter ao that appears in such analyses can be related to a 
moment of the coupling in the infrared. On the other hand our new parameter /i£ is more 
related to the 'reshuffling' of momenta associated with the production of mass — in this 
sense it is a much more intimately associated with the dynamics of hadronisation, and there 
is perhaps even a possibility that its determination could give us qualitative information 
about hadronisation. It is not currently clear whether perturbative contributions could also 
give contributions with a similar Q-dependence. 

Our analytical predictions for differences between different measurement schemes agree 
well with results from two commonly used Monte Carlo event generators, Herwig and Ariadne. 
Furthermore, in fits for as and ao from a range of event-shape variables, they explain a pattern 
observed both in e"'"e~ and DIS of the jet masses giving significantly larger values of both as 
and ao than other variables: this is because in the default schemes there are positive (non- 
universal) mass-corrections for the jet masses and negative corrections for the other variables. 
Measuring all variables in the ii^-scheme leads to a significant improvement in the consistency 
between the jet masses and the other variables. We also note that for the jet masses the p 
and E'-schemes are relatively insensitive to certain experimental systematics (associated with 
difficulties in identifying hadrons) which are relevant in the default measurement scheme. 

The Bad. Unfortunately it seems that a direct experimental verification of the properties 
of mass effects is quite difficult. In principle mass effects can be seen by looking at the differ- 



ence between two measurement schemes for the same variable (say massive and i?-schemes in 
the case of the jet mass). However the experimental determination of such differences relies 
on the accurate identification of the mass of each particle in an event, whereas the experi- 
mental procedure usually just involves the assignment of the pion mass to each hadron. Since 
a large part of mass effects seems to come from kaons and nucleons this is a poor approxima- 
tion. Given that the results are then corrected for what cannot be seen, using Monte Carlo 
programs, the resulting 'measurement' of the difference between two schemes is likely to be 
as much a reflection of the properties of the Monte Carlo events as of the actual events!|3 

It so happens that in one of the Monte Carlo programs most commonly used for correcting 
data, Pythia, the Q-dependence of mass effects differs significantly from our predictions 
at very high energies, going very roughly as (lnlnQ)/Q rather than {hiQ)'^/Q. This is 
probably related to approximations in Pythia's implementation of coherence. Fortunately for 
measurements in the p or ii^-schemes this should not have too large an effect on the practical 
determination of event shape values since the 3-momentum or energy are close to what is 
truly measured by detectors (it is the difference between them that is poorly measured). 
Furthermore at today's energies the discrepancy in the Pythia energy dependence is a small 
effect compared to other differences between event generators. 

In general, given the difficulties of a direct measurement of mass effects in event shapes, 
for an experimental test of the picture outlined here it might be worth investigating the 
feasibility of measuring some other observable expected to depend on the same (—1)**^ moment 
of particle energies. The simplest might be ^ • 1/ Ei where the sum runs over charged tracks. 
This might give an idea of whether the LPHD hypothesis used in this paper works also for 
negative moments, where it has so far never been tested. 

One 'negative' implication of our results relates to the determination of the perturba- 
tive and non-perturbative parameters of QCD. We would like to measure our new non- 
perturbative parameter /i^. But the degeneracy in a 3-parameter fit for as, and /i£ is such 
that the currently available data are not precise enough to give us any meaningful constraints 
on any of the fit parameters. Furthermore the fit results are very unstable with respect to 
systematic uncertainties. Given this limitation we might decide that for the time being we 
should carry on as before with two-parameter fits and ao (of course with the event-shape 
variables now measured in the £^-scheme). Our having neglected a parameter in the fit will 
translate to systematic errors on and a^. We can try to gauge the size of the systematics 
by carrying out the fit on results at different levels of the hadronic decay chain (using Monte 
Carlo results to determine the correction) — this suggests systematics of the order of ±0.004 
on ctg and of about ±0.1 on uq. 

The Ugly. In the course of our studies we have come across two effects, unrelated to light- 
hadron masses, that have a significant impact on event-shape studies and so warrant further 
investigation. 

Firstly there is the observation that heavy-quark effects have a large impact on fit values 
for Us and ao — measuring an event shape before and after heavy quark decay modifies as 
values by up to 0.006, even though at Mz the effect of heavy quark decay is only 2%. This 
pecuhar behaviour is seen because the / contribution from quark decay is simulated by 
an increase of and a decrease of a^. Given that the tools for studying heavy-quark effects 
in event shapes are well-established, there is a strong argument for carrying out analyses that 
make full use of them. 

^°One possible more direct experimental test of the ideas at the base of the predictions for mass effects 
would be a measurement of a quantity such as the sum of the inverse energies of all nucleons (or some other 
particle species). This should scale as (InQ)^. 



Secondly fits for tlie lieavy-jet mass (a very non-inclusive variable) lead to values for 
which are about 10% smaller than for inclusive variables like the thrust or the mean jet 
mass. This needs to be understood. It could be due to a difference in the behaviour of the 
perturbation series at higher orders. But in appendix there is evidence from Monte Carlo 
simulations that hadronisation corrections for ph have unusual characteristics: in contrast to 
what is seen in more inclusive variables, the hadronisation depends strongly on the underlying 
hard configuration. There is therefore a need to develop techniques allowing a more formal 
approach to the study of such problems. 
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A Summary of notation 

For convenience we give here a summary of the definition of the various schemes introduced 
in this article. 

p-scheme Scheme in which the observable is defined solely in terms of particle 

3-momenta. 

i^-scheme Scheme in which the observable is defined solely in terms of particle 

energies and angles. 

decay-scheme Scheme in which all massive particles are decayed isotropically into 
pairs of massless particles. The observable is then calculated using 
the resulting ensemble of massless particles. 

We also summarise some of the other notation used and introduced in this article. 
V An event-shape variable. 

Vp, Vei Vdecay An evcut-shape variable in p, E or decay-scheme, respectively. 

Cv The coefficient of the 'traditional' power correction for the observable 

V, introduced in eq. ( |3.5| ) and given for a range of observables in 
table |. 

{^mV) The non-universal mass-dependent correction to the mean value of 

the observable V, cf. eq. p.lU] ). 

5c^{m? / k^) The modification to cy for a particle with a given m^/ kf, cf. eq. ( ^.121 ). 



7v 



The coefficient of the leading ni?/kl dependence of 6c\;{m'^/kf) for 
fcf ^ m, as defined in eq. ( p.l8| ). It is shown for a range of variables 
in their default schemes in table 0. 



Xdecay(^^/fct ) The mean relative change in the sum of \kt\^s due to the decay of a 
massive particle of mass and transverse momentum k^, as defined 
in eq. (|3^ . 

A Shorthand for ACa/Po- Mass- dependent corrections have a leading 

term proportional to {\iiiQ/A)^/Q. 

El The normalisation of the non-universal mass-dependent correction for 

a hadronic level c/. eqs. ( ^.271 , ^.281 ). 

ixi The normalisation of the universal mass-dependent correction for a 

hadronic level £, c/. eqs. (|4.28|, |4.29|). 



B Broadenings in the decay scheme 

In determining the QCD-based predictions for mass-effects in section [4.2| we made the as- 
sumption that for all variables considered, mass effects matter only in the region of small 
rapidities. As a result we could ignore the rapidity dependence of the hadron distribution, 
because it is significant only at large rapidities. More precisely we expanded the rapidity 
dependence of ^h-wijl-, ^) and showed that it gave terms suppressed by a power of Og, 



y d^bh^Mols^ = 0{a,), (B.l) 



which holds if Sfy^^jiv) decreases sufficiently rapidly with 77. This works for all variables 
except the broadenings in the decay scheme, the reason being that for the broadenings the 
difference between the normal and decay schemes is sensitive to particle masses for particles 
up to rapidities of order For the example the total jet broadening has 

'^/BT,dccay,^(^) = fsAv) ^dccay,^ (B.2) 

where we can write iTSl 



(We recall that here 7]q has been defined with an additional 3/4 compared to what is given 
in |]15| ) . Since fsj. is roughly 1 up to 77 ~ 1 / ^/a^ the integral ( |B.1|) is of order 1 (while the 
leading term is of order cf. table |l|). Even though this term is subleading, in analogy 

with what is done for the cy's, we wish to control it, which can be achieved as follows. First 
we write 

J d7] 6fv,M ^hMv, Y) = $„,-.(0, Y)6cv,^ (B.4) 
where the r] dependence of ^h,-uj has been absorbed into 6cv^^: 



Substituting in the above equations for ^/sydocay,'^ gives 



In section P]l we had stated that 

If we compare our result (|B.6|) with cbj. in table |l|, we find that there is an additional term 
~4^^ = — not present in cy. This is the piece which arises from the rapidity dependence 
of the hadron distribution and in practice it gives quite a large correction to 5cbt^^^._^ 
The corresponding derivation for the wide-jet broadening involves the use of 

^fSwA^.^yAv) = Ibw iv) -^decay.L^ , (B.7) 

where 

We thus obtain 

jr _y ( /30 ^0 ^ /^O \ ,-p,Q^ 

The final form for the broadening mass- dependent power correction is analogous to 
eq. ( [4.26| ) but with the appropriate cj-dependence introduced for 7v: 



with 



(5,n5dccay) = 7i?...ay,A ^ J + ' ' ' (^-10) 



and 



7Bw,dccay,- 4 l^4v/2C^ 24Cf 2 ^ 16Cf; ' ^ ' ^ 

C Heavy quark decay 

While effects due to heavy quarks are not strictly speaking the concern of this article, it turns 
out that they can have a significant effect on the fit results for and a^. The presence of 
heavy quarks affects both the perturbative and non-perturbative contributions to the event 
shapes. 

Perturbative calculations of event shapes involving heavy quarks have been in existence 
for a few years now p8|-|4^. They have started to be used for experimental studies of 
events with heavy quarks, with hadronisation corrections deduced from Monte Carlo event 
generators (see for example [^). 

Power corrections to event shapes with heavy quarks have been studied in []T^ — they 



differ from the light quark case because very collinear radiation {rj > InQ/m^, where niq is 
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Table 6: Corrections to event shapes arising from heavy quark decay. 

the heavy-quark mass) is suppressed. For event shapes hke r, C and the jet masses this leads 
to a reduction of the power correction by an amount of order m^A/Q^. The broadenings are 
more complex. 

For typical measurements of event shapes, in the presence of heavy quarks there is a 
second 'hadronisation' contribution due to the fact that what is measured is not the final state 
involving charm or bottom-quark hadrons, but rather a final state where the heavy-quark 
hadrons have decayed. This has been discussed in some detail in the context of fragmentation 
functions in P^, and many aspects turn out to be quite similar for event shapes. Since (for 



Q ^ rrig) the quark is fast-moving, the effect of the decay is to produce a bunch of nearly 
coUinear hadrons whose invariant mass is the heavy-hadron mass (for simplicity, from here 
on we neglect the distinction between the heavy-quark and the heavy-quark hadron, while 
bearing in mind that the rUg relevant for the decay is actually the hadron mass). 

If we consider a Born configuration consisting of two back-to-back heavy quarks then it 
is quite straightforward to see what effect the decay will have on the simpler event shapes. 
In the massive scheme p is unchanged by quark decay. In the p or i?-schemes it goes from 
being zero before the quark decay to 

2 

mf, 

after quark decay. This result can also be used to deduce the correction to the thrust and 
C-parameter because in the collinear limit we have |H6[ 



tc^Pl + Pr Cc^6{pl + Pr). (C.2) 

where p^ and pr are the left and right-hemisphere jet masses respectively. To summarise, the 
corrections expected as a result of the decay of primary heavy quarks are shown in table ^. 

For the jet broadenings the situation is more complex. If we start from the Born config- 
uration and let the heavy quarks decay, then with respect to the thrust axis all the decay 
products will have transverse momenta of order ruq, leading to values for the broadenings of 
order rriq/Q. However to consider just the effect of quark decay on the Born configurations is 
an oversimplification. The variables discussed above had the property that they are linear for 
soft and/or collinear particles (even ph, to within the approximations required here): namely 
the effect of quark decay is the same regardless of whether we have the Born configuration, 
or one with extra soft and/or collinear particles. 

But the broadenings do not have this property. What goes on with heavy-quark decay is 
very similar to the dynamics that led to the rather complex form for the power correction to 
the broadening. Essentially the extra transverse momentum from the decay only contributes 
at order rUq/Q if the quark's angle with respect to the thrust axis is less than niqlQ. If the 
quark's angle is larger than this, then azimuthal averaging causes the mq/Q correction to 
be reduced to a rri^/ptQ correction, where pt is the quark's transverse momentum. A more 
quantitative understanding would require a full treatment of the 3-body heavy hadron decay 
and a study (from a perturbative resummation) of how the quark transverse momentum 
compares with rriq as a function of Q (for very large Q it will usually be much bigger, but 
phenomenologically accessible values of Q may not be large enough). 
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Figure 11: The hadronisation correction as a function of the value of the variable at parton 
level. 



In it has been pointed out that for the longitudinal fragmentation functions there 
can also be corrections proportional to rub/Q associated with the decay of secondary heavy- 
quarks produced from a soft gluon, though it is suggested that for today's energies such a 
behaviour may not yet have set in. The possibility of a similar contribution in event shapes 
should be investigated. 



D The heavy- jet mass 



We observed in section ^ that even in a 'proper' scheme the (as, ao) fits for the heavy-jet 
mass (and perhaps also the wide-jet broadening) seem to some extent inconsistent with the 
results for the other variables. The distinguishing feature of the heavy-jet mass is its non- 
inclusiveness, since it measures a specific hemisphere of the event (the heavy one), whereas 
other variables measure the properties of the whole event. 

We may well ask why non-inclusiveness leads to differences. One interesting analysis 
has been presented in |T^, which suggests that hadronisation corrections can be different in 
the two hemispheres and convert a perturbatively light jet into a heavy one. However this 
effectively increases the power correction rather than decreasing it and so cannot explain 
the relatively small ag and ao values that are observed. This does not mean that such a 
mechanism is not present at all — indeed in the difference between the E and p schemes the 
heavy jet mass correction is larger than that for the single jet mass [of. fig. and this could 
be due to such a mechanism (it could also simply be because there are more hadrons in the 
heavy hemisphere). 

To help understand what is happening we have used Pythia to look at the mean hadro- 
nisation as a function of the value of the variable at parton level, figure ITT|. For the thrust, 



the hadronisation is fairly independent of the parton- level thrust value. For the heavy jet 



masses (in the massive and j9-schemes) however the hadronisation correction is very negative 
for larger values of ph at parton level (we note that there are also very large event-by-event 
fluctuations in the hadronisation — this means that the hadronisation does not just cause 
a simple shift of the perturbative distribution). This feature may be at the root of the 
non-universality seen in the heavy jet mass, and needs to be understood — one explanation 
might be that when one hemisphere is perturbatively heavy, the phase-space that remains 
for further emissions is limited, and the only way of emitting non-perturbative radiation in 
the light-hemisphere is if some energy (and mass) is removed from the heavy hemisphere. 
However such a hypothesis would need to be placed on a more rigorous mathematical footing 
for it to be tested. 

Such a phenomenon might also play a role in the wide-jet broadening, where it has been 



observed that the theoretically predicted distribution is too wide at larger values of Bw [47 
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